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Abstract 

D ' It is shown that the density-potential mapping and the V-representability problems in the time- 
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dependent current density functional theory (TDCDFT) are reduced to the solution of a certain 



many-body nonlinear Schrodinger equation (NLSE). The derived NLSE for TDCDFT adds a link 
which bridges the earlier NLSE-based formulations of the time-dependent deformation functional 
theory (TDDefFT) and the time-dependent density functional theory (TDDFT). We establish close 

O ' 

relations between the nonlinear many-body problems which control the existence of TDCDFT, 
TDDFT, and TDDefFT, and thus develop a unified point of view on the whole family of the 

> 

\Q 1 TDDFT-like theories. 
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I. INTRODUCTION 



An accurate description of quantum many-electron dynamics is one of the most challeng- 
ing problems in the modern theoretical physics. Because a direct application of the standard 
many-body techniques to interacting nonequilibrium systems becomes prohibitively compli- 
cated, the time-dependent density functional theory (TDDFT)[l| remains practically the 
only available tool for modelling dynamics of realistic systems at a reasonable computa- 
tional cost. By formally avoiding a full solution of the complicated many-body problem, 
TDDFT gives a direct access to experimentally relevant observables, such as the density of 
particles or the density of current. Within the most common Kohn-Sham (KS) implemen- 
tation of TDDFT the exact density is calculated by propagating mean-field-like equations 
of motion for a fictitious system of noninteracting KS particles, which enormously simplifies 
the problem. Therefore a growing popularity of TDDFT is absolutely not surprising (see, 
e. g., recent special issues of PCCP [2] and TEOCHEM {3)], as well as this volume). 

Practical applications of TDDFT to realistic many-electron systems unavoidably employ 
approximations. However, in almost every paper, where TDDFT is applied, it is introduced 
as a formally exact theory. Unfortunately, at the present stage, the common statements of 
exactness express our expectations/believes more than the actual situation in the field of 
mathematical foundations of TDDFT and related theories. At the formal level TDDFT is 
based on two main statements: (i) A map from the external potential to the density is unique 
and invertable; (ii) A given density can be obtained from the Schrodinger dynamics driven by 
a properly chosen external potential (for example, the density coming from the interacting 
system can be reproduced by adjusting a potential in a fictitious noninteracting system). 
The first statement is known as the Runge-Gross mapping theorem J4]. The statement (ii) 
is commonly referred to as a V-representability problem, which is a corner stone of the KS 
formalism. 

A disturbing fact is that the Runge-Gross proof |4j of the mapping theorem is only valid for 
potentials which are analytic functions of time (t-analytic) and can be represented around 
the initial time by a convergent Taylor series. The conditions of the V-representability 
theorem, as formulated by van Leeuwen, are even more restrictive. One also assumes the 
t-analyticity of the density to uniquely recover coefficients in the power series expansion of 
the potential supporting that density, while a uniform convergence of the recovered series 
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is taken for granted. Therefore our present justification of TDDFT and its generalizations, 
such as time-dependent current density functional theory (TDCDFT) , is limited to a 
very narrow class of external potentials. It should be emphasized that the t-analyticity of 
the potential, and even its infinite smoothness in space, are not sufficient to guarantee the 
i-analyticity of the density [8j required for the present formulation of the V-representability 
theorem [5[. Hence, strictly speaking, the precise conditions of V-representability are still 
unknown. Obviously, a better understanding of the fundamental issues of TDDFT is highly 
desirable, especially in a view of explosively growing number of practical applications. 

A commonly accepted weak point of the existing proofs both in TDDFT 0, Q and in 
TDCDFT Q] is that the mathematical statement of the problem explicitly relies on the 
power series representation of the densities and/or potentials. It is absolutely clear that 
no further progress is possible without finding a new formulation, which does not require, 
at least in principle, the power series expansion of the basic quantities. In several recent 



works 



12| such a new route to the proof of existence of TDDFT-like theories has been 



proposed. Namely, the density-potential mapping and the V-representability problems can 
be restated in the form of the existence and uniqueness of solutions to a certain universal 
nonlinear Schrodinger equation (NLSE). The universal NLSE has been first formulated in 

nn 

the context of the time-dependent deformation functional theory (TDDefFT) [9|, [10j, where 
it appears as a natural step in solving the many-body problem in the comoving reference 
frame. In Ref. |8[ a conceptually similar formulation has been proposed for TDDFT, and 
more recently the NLSE approach has been adapted to rigorously prove the uniqueness and 
existence theorems for a lattice version of TDCDFT [ill ]. 

Despite the NLSE formulations proposed for TDDefFT in jt], [h]] and for TDDFT in 
js] are similar conceptually, their formal realizations look very different and appear to be 
completely disconnected. One of the aims of the present paper is find a relation between 
these two theories. I will derive the NLSE formulation of TDCDFT and show that it is 
precisely the missing link which bridges the two abovementioned proposals. We will see that 
the concept of a universal many-body NLSE provides us with a unified point of view on the 
whole family of TDDFT-like formalisms, namely, TDCDFT, TDDefFT, and TDDFT itself. 
This unification is the main result of the present work. 

The structure of the paper is the following. In Sec. IIA I formulate the NLSE approach 
to the mapping and the V-representability problems in TDCDFT. An equivalence of the 
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derived NLSE for TDCDFT and the NLSE for TDDefFT Q, Q is proven in Sec. IIB. In 
Sec. IIC I establish a connection of the NLSE formalism to the Vignale's approach to the 
TDCDFT mapping theorem [zj]. In Sec. Ill I derive a new NLSE formulation of TDDFT, 
and demonstrate its connection, on one hand, to the TDCDFT-NLSE of Sec. II, and, on 
the other hand to TDDFT-NLSE proposed in Ref. [8]. This uncovers interrelations between 
all members of the TDDFT-family and completes the unification. The main results of the 
paper are summarized in Conclusion. 

II. TIME-DEPENDENT CURRENT DENSITY FUNCTIONAL THEORY 
A. Nonlinear inverse many-body problem in TDCDFT 

Let us consider a system of N identical particles in the presence of time dependent 
external scalar C/(x, t) and vector A(x, t) potentials. The corresponding many-body wave 
function ^(xj., . . . , Xjy, t) is a solution to the time-dependent Schrodinger equation 

i^*(xi, . . . , xjv, t) = HVfa, ...,x N ,t) (1) 

with the following Hamiltonian 
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where Vj = d/dxj, and V(|x — x'|) is the interaction potential. For a given initial condition, 

*(X!, . . . , Xjv , 0) = * (xi, • • • , x^), (3) 

the dynamics of the system is completely specified by Eq. ([I]). 

The Schrodinger equation ([1]) is invariant under the following gauge transformation 

tf(t) e-'S;^.')*^), (4) 
A(x,t) -> A(x,t) + Vx(x,t), (5) 
f/(x,t) ^ U{x,t) - dt X (x,t) (6) 

Because the physical results are independent of the gauge we are free to choose a gauge 
function x(x, t) which is most suitable for a particular problem. We will see that TDCDFT 
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is most elegantly formulated in a temporal gauge that corresponds to the following choice 

a 

X(x,t)= [ U(x,t')dt'. (7) 



The transformation with the gauge function (jTJ) eliminates the scalar potential so that the 
many-body Hamiltonian reduces to the form 

g [Al = f '-' V --^- i » 2 + i^ y( |x J -x t |). (8) 

j=l j^k 

The key physical observables in TDCDFT are the density of particles n(x, t) and the 
current density j(x, t): 

n(x,f) = (¥(f)|n(x)|¥(f)), (9) 

j(x,t) = (*(t)|f (x)|*(t)> - ^U(x,t), (10) 

m 

where is the solution to the Schrodinger equation (OQ), and n(x) and j p (x) are the 
operators of the density and of the paramagnetic current, respectively, 

N 

n(x) = ^(x-x,), (11) 
j'=i 

. AT 

j' p w = £E{ v ^ 5 ( x -^)}- ( 12 ) 

The gauge invariance implies that the number of particles is locally conserved, i. e. the 
density and the current are connected by the continuity equation 

Mx,() + Vj(x,t) = 0. (13) 

This, in particular, means that n(x, t) is uniquely determined by the current j(x, t) and the 
initial density distribution n(x, 0) = ^o( x )- 

The usual statement of the problem in quantum mechanics corresponds to solving the 
time-dependent Schrodinger equation flTJ), (jSJ), for a given initial condition ([3]) and the ex- 
ternal driving potential A(x, t). This determines the many-body wave function ^[^o A](i) 
and thus any observable, e. g. the current j[\&o> A](x, t), as unique functionals of the initial 
state \l/o and the vector potential A. The existence of TDCDFT assumes that the inverse 
problem also possess a unique solution. Namely, TDCDFT is valid if, given the initial state 
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\l/o and the current j(x, t), one can uniquely reconstruct the time- dependent wave function 
^[^ch j](t) and the potential Af^o, j]( x ; t) which supports the prescribed current. 

The key observation is that, by "reinterpreting" the definition fllOp of the current, the 
above inverse problem can be mathematically posed in a form of the following nonlinear 
system of equations 

id&(t) = H[A]V(t), (14) 

777 r 

A(x,t) = — — (*(*)|f(x)|*(*)>-j(x,t) , (15) 



which has to be solved with the initial condition ([3]). The Hamiltonian in ( TT4")) still has a 
form of Eq. ([8]). However, now the vector potential A is not fixed externally, but calculated 
selfconsistently from Eq. (1T51) . The density n(x, t) entering (|T5l) is determined either by 
intergating the continuity equation ({TBI or calculated directly as the expectation value of 
Eq. (Q. Both definitions are identical if the wave function \^(t) satisfies Eq. f|T4l) . By insert- 
ing the selfconsistent vector potential (|T5l) into (fT4j) we obtain NLSE with a special cubic 
nonlinearity. The solution of the Cauchy problem for this NLSE, provided it is well posed, 
returns the time-dependent many-body wave function and the vector potential as universal 
functionals of the initial state \l/o and the given current j(x, t). Therefore the problem of 
existence of TDCDFT reduces to proving the existence and uniqueness of solutions to NLSE 
defined by Eqs. and ffi~5]) . 

To illustrate how the NLSE approach works I consider the simplest case of one quantum 
particle. For N — 1 the nonlinear Cauchy problem (JL41 . (|T5l) . and ([3]) simplifies as follows 

i$tt(x, t) = 2^HV - A(x, t)) 2 M/(x, t), (16) 



A(x,t) 



1^1 



_!(^*V^ - W^*) - mj(x, t) 



(17) 



*(x,0) = vp (x) = V^W^ (X) . (18) 

On the first sight this nonlinear problem may look complicated, but it turns out to be 
trivially integrable. The corresponding exact analytic solution can be represented in the 
form 



tt[*o,j](x,*) = vMx,t)e^\ (19) 
A[* J](x,t) = V^(x,t)-m^4, (20) 

7i(X, t) 



where the functions n(x, t) and y?(x, t) are defined as follows 
n(x,f) = no(x) - / dt'Vj(x,t'), 
p(x,t) = V9 (x) + 



dt' 



V 2 v^(xTtO mj 2 (x,t') 



2my / n(x, £') 



2n 2 (x,t' 



(21) 
(22) 



Equations ( H~9]) — f[22T) provide us with an explicit example of the universal TDCDFT func- 
tional ^[^oj j](t) and Af^o? j]( x , t) recovered from the NLSE (fl4j) . (fl5j) for a particular case 
of N = 1. 

Apparently the general NLSE for any iV > 1 is not solvable analytically. However the 
mathematical structure of equations remains practically the same, which leaves a strong 
hope that early or later the uniqueness and existence theorems for the general nonlinear 
many-body problem, and thus the existence of TDCDFT, will be proved. In fact, recently 



I have succeeded to prove the corresponding theorems for a lattice TDCDFT 



ll| which is 



formulated in terms of a discrete (actually finite-difference) version of the above NLSE. 



B. Equivalence of the nonlinear universal problems in TDCDFT and TDDefFT 

The NLSE formulation has been first formulated in the context of TDDefFT where it 
appears as a natural "universal" step in solving the many-body in the comoving reference 
frame 0, [h]] . Below I will show that there is a one-to-one correspondence between the uni- 
versal problem derived in Refs. @,Q for TDDefFT and the NLSE (fTlh (1150 for TDCDFT. 
To prove this correspondence we proceed as follows. First, for a given current j(x, t), we 
integrate the continuity equation ( jl~3]) to calculate the density n(x, t), and then define the 
velocity field: 

= < 23 > 

Using the velocity field v(x, t) one finds a set of Lagrangian trajectories x(£,t) by solving 



the following initial value problem (see, e. g. Refs. [9l. llOj) 



*(£,<) = v(x(£,f),*) x(£,0)=£, (24) 

where x(£,t) = <9jx(£,£). Physically the function x(£,i) is a trajectory of an infinitesimal 
fluid element which start its motion at t = from the point £. The formal significance 
of the trajectory function x(£, t) is that the equation x = x(£, t) defines a transformation 



of coordinates x — y £ which corresponds the transformation from the laboratory frame 
(x-space) to the comoving Lagrangian frame (£-space). 

The last step in proving the equivalence of NLSE in TDCDFT and TDDefFT is to 
rewrite Eqs. (|14p and f lT5"j) in the comoving frame. Formally we perform the transformation 
of coordinates x — y £ : x = x(£, t) and define the transformed many-body wave function 



^(£i> ■ ■ ■ ) £n, t) in the comoving frame as follows (I], 



3 



N 



...,tN,t) = Y[g* t)e- iS ^^{^ *),..., x(£ JV , t), t), (25) 

3=1 

where g(£,t) = det[g Mi ,(£, t)] is the determinant of the metric tensor induced by the trans- 
formation of coordinates (y/g is the Jacobian of the transformation) 

9r&t) = ^*) = W" 1 = ^if ( 26 ) 

and S c i(£,t) is the classical action of a particle moving along the trajectory x(£, t) 

5d(t = jT [y (*& i)) 2 + t)A(x(& . (27) 
Note that the factor FIjli in PS) 

ensures the standard normalization of the wave 
function = 1 after a non- volume-preserving transformation of coordinates. After 

straightforward algebra we find that the transformed Schrodinger equation ( I14p and the 
"self-consistency" equation ( FlBTl can be reduced to the following form 

...,£ N ,t) = H[ 9ij , ...,£ N ,t), (28) 

777 a ■ — - 

A^,*) = —j-Amimmt))- w 

"o(«) 

In Eq. (128]) H[gij,A\ is the Hamiltonian (jSJ) transformed to the comoving frame: 

H[ 9lJ ,A] = ^^J^K^ + (30) 
i=i ~~ Mi 

where A^ iM = — ic^ — A^($j,t), g^ v = g^ u (£j,t), and l^. is the distance between jth and 
fcth particles in the moving frame (the length of geodesic connecting points £j and in the 
space with metric g^)- The selfconsistent vector potential A.(£,t) in ff28l . fl29|) is related to 
A(x, i) entering (TT^j) and (|T5|) as follows 

MM = Q^Mx(U),t) + ^ - fys^t). (si) 
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The first term in the right hand side in this equation is the transformed left hand side in 
f|T5|) . The second term originates from the term mj/n in f|T5|) . while the last term in ( l3Tj) 
comes from the expectation value of the paramagnetic current in and is related to the 
extra phase factor in the transformed wave function f l2"5"j) . 

The nonlinear system of equations (I2"8"j) . ( T2§|) is exactly the universal many-body problem 
formulated within TDDefFT (see Eqs. (58), (59) in |a| or Eqs. (3.38), (3.39) in 10|). In 
fact, the NLSE in TDCDFT and TDDefFT correspond to the same mathematical problem 
formulated in the laboratory or the comoving frame respectively. Hence to prove the exis- 
tence of TDCDFT we are free to choose any of the two formulations. The advantage of the 
comoving frame formulation is a much simpler form of the selfconsistency equation. Indeed 
Eq. ( 12"§|) does not contain an inhomogeneous term and involves only the initial density dis- 
tribution uq. The price for that - the appearance of the time-dependent metric tensor in 
the Schrodinger equation - is probably not too high from the conceptual point of view. On 
the other hand, the formulation derived in this paper, Eqs. ( fl4l) and ( !T5l) . is much easier 
to digest as it is based on the usual quantum dynamics in the laboratory frame. 

C. Connection of the nonlinear inverse problem in TDCDFT to the Vignale proof 
of the mapping theorem 

In Ref. {t| Vignale generalized the van Leeuwen construction Q to prove the mapping 
theorem for TDCDFT. The main idea of the proof is to relate the vector potential A(x, t) 
to the current and the many-body wave function using the force balance equation. Then 
this relation has been interpreted as an equation of motion for A(x, t) and solved recursively 
assuming t-analyticity of the vector potential. The present subsection is aimed at demon- 



strating that the statement the problem used in [7] directly follows from the nonlinear system 
of Eqs. (HU) and (fT5|). 

Let us first multiply both sides of (|T5|) with the density n(x, t) and differentiate the result 
with respect to time: 

nd t A + Ad t n = md t (*(t)|j p (x)|*(t)) - md t y (32) 

The time derivative of the paramagnetic current j p (x, t) = (\l/(£)|j p (x)|\l/(i)) in the right 
hand side of ( )32l) can be straightforwardly calculated using the Schrodinger equation (fl4l) . 
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The result of the differentiation takes the following form 

mdtj* = A^dtn + \(r - -A) x Bl - d u U^ u , (33) 

L V Tfl / J fj, 

where B = V x A is the magnetic field associated with the selfconsistent vector potential, 
and the last term is the stress force with n^(x, t) = (\I/(t)|fl^(x)|\l/(t)) being the stress 
tensor. An explicit form of the stress force and the stress tensor can be found, for example, 
in 0, H, 13]. Inserting the result of ( 133]) into ( 132]) we arrive at the following equation of 
motion for the vector potential 

ndtAp - 



Up _ ™ A ) x Bl - ctn M „ - md t3 „ (34) 

which is exactly the force balance equation. Using the definition of ffTO]) . which is an ab- 
solutely legitimate operation at the solution point, we find that ( l34j) becomes identical to 
Eq. (10) in Hence the force balance equation, which has been used in [3] to construct a 
recursion for A[j](i), is nothing but the time derivative of the selfconsistency equation ffl~5]) . 
In other words, the force balance equation, considered as the equation of motion for the 
vector potential A(x, t), can be integrated explicitly and its solution is given by (fl~5]) . 

Apparently the equation for A ( fl~5l) is advantageous as it leads to a much simpler and 
cleaner form of the nonlinearity in the inverse many-body problem. Another important 
point is that Eq. ( Ti~5l) is local in time. Just from the structure of the nonlinear problem (Tl4l) . 
(|T5|) we immediately see that, if the solution exists, the vector potential A[j](t) is a strictly 
retarded functional of the current j(x, t). In particular, it can not contain even local terms 
which depend on the time derivative of the current. The time locality of (fT5|) also makes it 
absolutely obvious that the power series expansion of A(t) is unique, which is the essence of 
the Vignale theorem Of course, to construct such an expansion we have to assume, as 
usual, the t-analyticity of the current. I am not presenting here the explicit construction of 
the recursion because it is practically trivial, and basically reproduces the recursion for the 
selfconsistent potential «4.(£) in the universal problem (|28j) . (129]) for TDDefFT 0,[l(3|. Hence 
the present derivation of the nonlinear inverse problem (jUj) . f[T5]) can be also viewed as a 
considerable simplification of the proof for the Vignale theorem. Clearly the convergence 
of the series remains unproved as in all power series based approaches to the problem of 

finn 

TDDFT and TDCDFT [5|, 13, 121] . The hope is that a cleanly stated nonlinear initial value 
problem, ([14]) , (fl5]) . can be treated using more advanced mathematical tools as it has been 



done recently for a lattice version of the theory 
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III. TIME-DEPENDENT DENSITY FUNCTIONAL THEORY 



In Sec. II we have found that the problem of existence of TDCDFT is equivalent to 
that of existence and uniqueness of a solution to a certain time-dependent NLSE. We have 
proved a one-to-one correspondence of this NLSE to the universal nonlinear problem derived 
earlier for TDDefFT, and demonstrated its close relation to the Vignale's power-series proof 
of the TDCDFT mapping theorem. Recently a conceptually similar NLSE setup has been 
also proposed for TDDFT [8[]. In this section I will derive an alternative formulation of the 
TDDFT problem, and establish a connection between the results of Sec. II (i. e., the NLSE 
formulation of TDDefFT and TDCDFT) to the NLSE proposed in Ref. [8|. 

The main statement of TDDFT is that the density n(x, t) completely determines the 
many-body dynamics, provided the latter is driven by a scalar potential. Hence our starting 
point is the many-body Schrodinger equation ([1]) with a Hamiltonian that contains only the 
external scalar potential ClS db driving force: 



N r v 2 11 

H = H ~^ + U ^^ +2£n|x,-x fc |). (35) 



j = l L 

At the formal level this standard form of the Hamiltonian assumes that a given, purely 
longitudinal external field is described using a Coulomb gauge (VA = 0). 

Following the ideas of the previous section we perform a gauge transformation @ with 
the gauge function x( x )^) defined by Eq. ([7]). In other words, we switch from the origi- 
nal Coulomb gauge to a temporal gauge, where the scalar potential is absent, but instead 
the dynamics is driven by a purely longitudinal vector potential A(x, t) = Vx(x, t). The 
Hamiltonian in the temporal gauge takes the form 

m = E HVj ~ l^' t)? + \ E nfe - „d. m 

3=1 " j^k 

Obviously, the driving force is still determined by a single scalar function x(x, t) which 
equals to the potential Z7(x, t) integrated over the time, Eq. ([7j). 

Solution of the Schrodinger with the Hamiltonian H[x] yields the many-body wave func- 
tion \l/[\l/o, x]{t) which can be used to calculate any observable, for example the density or 
the current, as a functional of the initial state \l/o an d the driving "potential" x( x >^0- It is 
worth noting that the density n(x, t) is always given by the gauge independent Eq. (jH]), while 
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the formal expression for the current j(x, t) in the temporal gauge acquires a diamagnetic 
contribution: 

j(x,t) = (¥(f)|ftx)|¥(f)> - r^vx(x,f) (37) 

For the discussion below the following simple observation is of primary importance. By 
taking the divergence of Eq. ( 1371) and using the continuity equation ( TT3l) we find that the time 
derivative of the density, d t n(x,i), can be calculated directly as the following expectation 
value 

9tn(x,t) = - W)|Vj p (x)W)> + ^V[n(x,t)Vx(x,t)]. (38) 
The relevance of this equation for TDDFT should be obvious from the familiar looking 

n 

Sturm-Liouville operator |5[ in the right hand side (the second term). 

The validity of TDDFT means that, given the initial state \l/o and the density n(x, t), 
we can uniquely recover the time-dependent many-body wave function \I/ [\I/o ? n }{t) an d the 
corresponding potential which gives rise to the prescribed density. In the temporal gauge this 
inverse problem can be formulated in terms of the following nonlinear system of equations 

id&(t) = H[xMt), (39) 
V[n(x,t)Vx(x,t)] = m(^(t)\W(x.)\y(t)),+md t n(x,t). (40) 

where Eq. (J4"D]) corresponds to Eq. fl38l) reinterpreted as an equation for x(x, t). The system 
(|39|) . (140]) has to be supplemented with the initial condition \I/(0) = and solved for a 
given n(x, t). 

The Hamiltonian in the Schrodinger equation ( !39l is the same as in the usual linear quan- 
tum problem, i. e. it is defined by Eq. ( 1361) . However, the "potential" function x(x, t) in H[x] 
is determined selfconsistently from Eq. (HU1) for a given density. Mathematically the system 
of (|39|) . (140]) is equivalent to NLSE with a spatially nonlocal cubic nonlinearity. Provided the 
Cauchy problem for this NLSE is well posed, its solution returns the functionals \l/[\l/o; n ](t) 
and n ]( x ! t). A functional, which corresponds to the scalar potential in the original 

gauge, is calculated as U[^q, n](x, t) = d t x[^o, ^]( x , t) according to the definition (J7J) of the 
gauge function. Hence the proof of the density-potential mapping and the V-representability 
reduces to proving the uniqueness and the existence of solutions for the nonlinear problem 
(|39|) . (j40|) . In this setting the existence corresponds to the V-representability, while the 
uniqueness of a solution is equivalent to the Runge-Gross mapping theorem. 
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A connection of the above NLSE for TDDFT to the corresponding problem in TDCDFT 
(see Sec. II) is practically obvious. Indeed Eq. (140]) is nothing but a divergence of Eq. ffl5|) 
when the latter is restricted to longitudinal vector potentials of the form A = V%- Hence 
TDDFT can be considered as a particular case of TDCDFT for systems driven by a purely 
longitudinal vector potential. Since in this case the vector potential is uniquely determined 
by a single scalar function x(x, t), a knowledge of a single scalar collective variable, the 
density n(x, t), is sufficient to uniquely recover both the full many-body wave function and 
the driving potential. 

Let us now demonstrate the equivalence of Eqs. (1311 . fHDl) to the NLSE formulation 
proposed in 8|. The argumentation below is very similar to that used in Sec. II to prove the 
equivalence of our NLSE formulation ( 1T4"]) . ( TTH]) . and the Vignale's approach to the TDCDFT 
mapping problem. As a first step we differentiate Eq. ( 140"]) with respect to t: 



V[nVx] + V[nVx] = md t (y(t)\Vf(x)\y(t)) + mh. (41) 

Next, we use the Schrodinger equation ([3"9~]) to compute the time derivative of the quantity 
(\I/(t)|Vj ?3 (x)|\E'(t)) entering the right hand side in (14T]) . This is equivalent to simply setting 
A = Vx in Eq. (133]) and taking its divergence, which yields the following result 

m$(#(t)|Vj p (x)|#(t)) = V[hV X ] ~ d^U^. (42) 

Inserting (142"]) into (14 ip we obtain the following form of the selfconsistency equation, 

V[n(x, t) W(x, t)} = mn(x,t) - d^U^t), (43) 

where t/(x, t) = 9 4 x(x, t) [see Eq. ([7])] is the usual scalar potential. Apparently Eq. ( )43|) 
is equivalent to the equation (5) in Ref. [8J. More precisely, after the transformation to 
the Coulomb gauge, equations (1431) and (1311 identically reproduce equations (5) and (6) of 
Ref. jjj. This proves the equivalence of the present formulation of TDDFT problem and 
NLSE proposed in 



Now we can answer the question raised by the authors o 



between the particular NLSE derived for TDDefFT in 
TDDFT in 



Ref. 



,8 



- what is the relation 



and the one proposed for 



8j. Starting from the TDDefFT universal problem, (T2"g|) and (T2"5|) . we succes- 
sively perform the following transformations, (i) First, we transform the TDDefFT-NLSE 
of Eqs. (|28|) . (129]) from the comoving to the laboratory frame. The result is the TDCDFT- 
NLSE, ([14]) . ( rT5]) . (ii) Then, in the TDCDFT framework we consider only longitudinal 

13 



vector potentials by setting A = Vx in (fl4l) . (fl5]) . and taking the divergence of Eq. (fl5l) . 
This brings us to the TDDFT-NLSE in the form ([39]), gO]). (iii) Next, we differentiate flU} 
with respect to time and, as a result, arrive at the TDDFT-NLSE of ( I39p . (I43p . (iv) Finally, 
we switch from the temporal to the Coulomb gauge, and transform Eqs. f[3"9~]) and (|43p to 
the TDDFT-NLSE proposed in Ref. |SJ. 

IV. CONCLUSION 

The main result of the present work is a unified approach to the density-potential mapping 
and the V-representability problems in TDDFT, TDCDFT, and TDDefFT. The existence 
of each theory in the above family is mapped to the uniqueness and existence of solutions 
to a certain nonlinear many-body problem. We have uncovered close relations between 
the nonlinear problems for all three theories, and established links to all previously known 
approaches to TDDFT and TDCDFT. The analysis performed in this work shows a spe- 
cial role of TDCDFT in the whole problematics of the TDDFT-like formalisms. Indeed, 
starting the nonlinear universal problem, ( f!4|) . ( fl5|) . for TDCDFT one can easily generate 
the corresponding nonlinear problems for the other members of the family. The NLSE for- 
mulation of TDCDFT demonstrates the main features of the new approach in the most 
transparent form. It is free of the additional Sturm-Liouville-type nonlocality of TDDFT, 
and does not involve complications related to the time- dependent metrics as TDDefFT. This 
clearly indicate that first it makes a sense to concentrate on the existence and uniqueness 
theorems for the TDCDFT nonlinear problem of Eqs. ( [14"]) . f fT5|) . Provided this problem is 
solved the corresponding results for the other members of the family should follow almost 
automatically. 
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